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PACS 03 . 75 . Lm - Tunneling in Bose-Einstein condensation 

Abstract. - We propose a new scheme for observing Josephson oscillations and macroscopic 
quantum self-trapping in a toroidally confined Bose-Einstein condensate: a dipolar self-induced 
Josephson junction. Polarizing the atoms perpendicularly to the trap symmetry axis, an effective 
ring-shaped, double-well potential is achieved which is induced by the dipolar interaction. By 
numerically solving the three-dimensional time-dependent Gross-Pitaevskii equation we show that 
coherent tunneling phenomena such as Josephson oscillations and quantum self-trapping can take 
place. The dynamics in the self-induced junction can be qualitatively described by a two-mode 
model taking into account both s-wave and dipolar interactions. 
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Introduction. — Josephson effects are a signature of 
quantum coherence in macroscopic many-body systems. 
Firstly predicted and observed when two superconduc- 
tors were connected through a weak link [I], Josephson 
effects have also been experimentally observed in a va- 
riety of systems: in superfluid helium flowing through a 
sub-micrometer aperture [5] and through an array of nano- 
apertures [3111]; in contact-interacting Bose-Einstein con- 
densates confined in a double- well trap [SUB] and in an 
optical lattice [7J; and very recently in exciton-polariton 
systems in semiconductors [8] . All these systems are real- 
izations of a Josephson junction (JJ). Internal Josephson 
dynamics has been also experimentally observed between 
different hypcrfinc states of a spinor Bose-Einstein conden- 
sate in Ref. [§]■ In Bose-Einstein condensates (BECs), due 
to the nonlinearity introduced by the s-wave contact in- 
teraction, there appears a new phenomenon called macro- 
scopic quantum self-trapping [5J[TUJ[TT], characterized by 
the locking of most of the atoms in one of the two wells. 

Dipolar Bose-Einstein condensates (dBECs) [T21[T5] of- 
fer a new playground for the study of Josephson effects due 
to the anisotropic and long-range character of their inter- 
action |14| . Here we present a novel scenario for inves- 
tigating the coherence properties of quantum transport: 
a dipolar self- induced Josephson junction (SI J J). Such a 
junction is not directly created by an external potential, 
but rather based on the appearance of an effective ring- 
shaped double- well potential in a toroidally confined dipo- 



lar condensate due to the anisotropy of the interaction. 
Josephson dynamics in dipolar condensates [15j[T6] and in 
spinor dipolar condensates [17] has been addressed in the 
literature for external double well potentials, but this is 
the first time to our knowledge that a SIJJ is investigated. 

System description. — We consider TV = 5 x 10 4 
atoms of 52 Cr with magnetic dipolc moment fi — 6 [Ib (be- 
ing fiB the Bohr magneton) confined in a pancake-shaped 
toroidal trap [18] 

V t (r) = ^-(oj ± rl + u z z 2 ) + V cxp(-2ri/ a 2 ) , (1) 



where m is the atomic mass, r± = \J x 2 + y 2 is the dis- 
tance to the trap symmetry axis, and w±_ = 8.4 x 2n s _1 
and lu z = 92.5 x 2ir s _1 are the radial and axial har- 
monic trap frequencies. Vq = 30 huj± and tro = 2.1 a± 
are the strength and waist of the Gaussian beam that 
creates a hole in the condensate along the z axis, with 
a± = y/h/(mu}±) the radial oscillator length. We consider 
the dipoles oriented perpendicularly to the trap axis (say y 
axis), and a small and repulsive s-wave scattering length, 
a = 14 ob (with ob the Bohr radius), to ensure that the 
dBEC is stable but dominated by the dipolar interaction. 
In this configuration the axial symmetry imposed by the 
external confinement is broken, and the anisotropic char- 
acter of the dipolar interaction is enhanced [T5]. This 
effect has been also addressed in dipolar Fermi gases con- 
fined in very narrow rings |19] , and recently, in bosonic and 
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fermionic dipoles in a quasi-one-dimensional ring trap [3D] . 

The combination of the external toroidal trapping and 
the mean- field dipolar potential, Vd (r) , defines an effective 
potential, V c s(r) = Vt(r) + Vd(r), which has the shape 
of a ring-shaped double well. This structure is shown 
in fig. [1] (a) for the z = plane, with the two potential 
wells in the direction perpendicular both to the trap sym- 
metry axis and the polarization axis. Since the centered 
Gaussian potential introduces a strong repulsive barrier at 
r±_ = that prevents the atoms from tunneling through 
it, the double well structure arises in the azimuthal direc- 
tion. Figure [T] (c) shows the minimum effective potential 
along the azimuthal coordinate, <p, with the two minima 
at ip = and ip = tt and the two barriers at f = tt/2 
and if = 3ir/2. Note that it is thus not defined at a fixed 
r±, but r± slightly changes with ip, being however close 
to r_i_ ~ 2.5 a_i_. According to this effective ring-shaped 
double well, the atoms localize mostly in the attractive 
regions inside the wells, producing an azimuthal density 
dependence. This is seen in panels (b) and (d) of fig. [T] 
which show the density and the maximum density along 
if on the z = plane, respectively. 
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Fig. 1: Effective potential (a) and corresponding ground state 
density (b) in the z — plane. Minimum effective potential (c) 
and maximum density (d) in the z = plane as a function of 
the azimuthal angle p. 

The system resembles two dipolar condensates in a ring- 
shaped geometry that are coupled via two links. If the 
barrier height, AV e g, is large compared to the chemical 
potential, fx, the system behaves as two weakly linked 
condensates. The weak link condition can be reached by 
tuning the s-wave scattering length to small values, which 
should still be large enough to prevent spatial spontaneous 



symmetry breaking [T5]. We find that the ground state 
fulfills AV c g/jl = 1.1, which is at the limit of the weak 
link condition. The system can thus be thought of as a 
self-induced Joscphson junction. In the remaining of this 
letter we will show that it can indeed behave as a bosonic 
J J, presenting Josephson as well as self-trapping dynam- 
ics. Note here that this junction really consists of two 
coupled SIJJs, which in this scenario behave in phase in 
much the same way as the array of nano-apertures does in 
the experiments with helium [3]|4]. 

For the sake of completeness we will briefly describe the 
system for scattering lengths below a = 14 ag. For val- 
ues still close to this one, the effect of a reduction of the 
scattering length is translated into an increase of the ra- 
tio AV c ff/p,, but with the effective potential still having 
the shape of a double well. However, reducing further the 
scattering length to about a = 12 as a spatial symmetry 
breaking occurs. All the atoms occupy only one of the 
wells, and the system resembles an origin-displaced cigar- 
shaped condensate |18| . Symmetry breaking phenomena 
in atomic condensates confined in double-well traps have 
been widely addressed in the literature, both for contact 
interacting BECs [2TJI22] and dBECs H3HB]. They are 
thought to be a signature of a quantum phase transition, 
during which the state of the system is strongly corre- 
lated and entangled, forming a macroscopic Schrddinger- 
cat state [2"2"H2~41 r For dipolar condensates, such strongly 
dclocalized states have also been theoretically proposed in 
a three- well configuration |25| . 

Dynamics in the SIJJ. — We study the dynamics of 
this dipolar SIJJ within the mean-field framework by solv- 
ing the full-3D time-dependent Gross-Pitaevskii (TDGP) 
equation: 



2m 



V 2 +V t (r)+ g\V(r,t)\ 2 + V d (r,t) 







(2) 



where g = 4Trh 2 a/m is the coupling constant of the con- 
tact interaction. Note that the mean field dipolar poten- 
tial, Vd(r,i), depends on time through the wave function 
^(r, t). At each time step, it can be written in terms of 
the microscopic dipole-dipole interaction as 



V d (T,t) 
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(3) 



where r and r' are the positions of two aligned dipoles and 
9 is the angle between their relative position and the po- 
larization axis. The TDGP equation with the dipolar term 
becomes an intcgro-differcntial equation. To solve the time 
evolution we have used a Hammings algorithm (predictor, 
corrector, modifier) initialized by a fourth-order Runge- 
Kutta method. The dipolar term has been treated using 
Fourier transform techniques (see Ref. [25] and references 
therein). The spatial grid used in the computation con- 
tains 128 x 128 x 64 points with a spacing of 0.12, 0.12 and 
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0.08 a_i_, respectively for the directions x, y and z. The 
time step used is 1.6 x 10 _4 o;J 1 . 

Josephson effects in BECs are characterized in terms of 
two conjugate variables: the population imbalance and 
the phase difference between the two wells [Sl ITUlfTTj . 
From fig. [T] we define the left and right wells as the 
regions where x < and x > 0, respectively. The 
population imbalance and phase difference are given by 
Z(t) = (N L (t) - N R (t))/N and 0(t) = </> R {t) - (f) L (t), re- 
spectively, where A?£(m(t) corresponds to the number of 
atoms on the left (right) well, and 4>L(R)(t) the phase of 
the dipolar condensate averaged on the left (right) well. 
To prepare the system with an initial population imbal- 
ance Z(0) 0, we solve the stationary Gross-Pitaevskii 
equation in imaginary time |18j with a small tilting in the 
external potential, ensuring thus the desired asymmetry 
in the population of the two wells. We then let the system 
evolve in time without any tilting, according to Eq. ([2]). 

For small values of the initial imbalance the system ex- 
hibits Josephson oscillations between the two self-induced 
wells. Figure [5] shows the dynamic evolution of the pop- 
ulation imbalance and the phase difference corresponding 
to initial conditions Z(0) = 0.1 and (f)(0) = (solid lines). 
Z(t) and 4>(t) present sinusoidal oscillations shifted by 7r/2, 
and the time average of the population imbalance is zero. 
In this regime, the atoms tunnel periodically from the left 
to the right well and back. The wave function of the dBEC 
at each well remains coherent during this process, that is 
with a uniform phase. This is translated into a phase dif- 
ference that oscillates in time at the same frequency as the 
imbalance. 
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Fig. 2: Time evolution of the population imbalance (top panel) 
and phase difference (bottom panel) with initial conditions 
Z(0) = 0.1 and (f>(0) = 0. Solid lines are the TDGP results, 
and dashed lines correspond to the TMM. 

For a large initial population imbalance, the SIJJ enters 
the regime of self-trapping oscillations. This can be clearly 
seen in fig. [3l where the imbalance and the phase differ- 
ence are plotted as a function of time for initial conditions 
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Fig. 3: Time evolution of the population imbalance (top panel) 
and phase difference (bottom panel) with initial conditions 
Z(0) = 0.65 and 0(0) = 0. Solid lines are the TDGP results, 
and dashed lines correspond to the TMM. 

Z(0) = 0.65 and cf>(0) = (solid lines). In this situation, 
the time average of the imbalance remains close to 0.5 and 
the phase difference is unbounded (running phase mode). 
Although in this regime the atoms remain locked in one 
of the wells, there is still some tunneling of particles at 
a frequency higher than in the Josephson regime. These 
imbalance oscillations are non-sinusoidal but have a more 
involved structure, which hints at a rich self-trapping dy- 
namics. 

In external double well potentials, the transition be- 
tween Josephson and self-trapping regimes is related to the 
presence of a saddle point of the Hamiltonian in the (Z, (f>) 
space that separates closed (Josephson) and open (self- 
trapping) trajectories with constant energy [TT1 12"51 12"§] . 
Extending the analogy between external and self-induced 
double-well potentials, one can expect also in the present 
case a similar critical behavior separating the two dy- 
namical regimes. By solving the time-dependent Gross- 
Pitaevskii equation, besides the results shown in figs. [2] 
and [3J we have obtained self-trapping dynamics for the 
initial conditions Z(0) = 0.5 and </>(0) = 0, which is com- 
patible with the critical imbalance obtained with a two 
mode model, Z c = 0.43 when </>(0) = (see below). 

Figures [2] and [3] show that a toroidally confined dBEC 
behaves as a SIJJ, and that such a structure can sustain 
both Josephson and self-trapping dynamics. In this sense, 
the SIJJ is robust, since the dipolar interaction is strong 
enough to keep the double-well shape of the effective po- 
tential in time, even when the initial imbalance is large. 
We have therefore shown that a new class of systems exists 
where the double-well potential structure is self-induced 
by an anisotropic interaction (in this case the dipolar in- 
teraction), and that for these systems Josephson and self- 
trapping oscillations are predicted to occur depending on 
the initial population imbalance between the two wells. 
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Two-mode model. — To gain insight into the tunnel- 
ing dynamics obtained by evolving the TDGP equation, 
we have performed a two-mode analysis of the SIJJ, tak- 
ing into account both s-wave and dipolar interactions. The 
two-mode model (TMM) lies on the assumption that the 
dynamical behavior of a JJ can be fully captured by ana- 
lyzing the coherent dynamics between two spatially local- 
ized modes: the left and right modes, respectively 
and $fl(r). In this approximation, the condensate order 
parameter is written using the ansatz |10 | lll j 



(4) 



(r, t) = Vi(*)$i(r) + ^«(t)$fl(r) 



with = Sij, and the coefficients fulfill ipj(t) = 

y/Njfye***®, i,j = L,R. Note that in this ansatz the 
time evolution is contained only in the coefficients tpj(t). 

By substituting ansatz (j4|) into eq. ([2]) and performing 
some algebra retaining all the overlaping terms [27f28] , one 
obtains the two-mode equations for a symmetric dipolar 
SIJJ, 



Z =(-l + a)Vl - Z 2 sin< 
Z 



--AZ-(-l 



with 
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(6) 

(7) 
(8) 
(9) 



The different parameters appearing in (O-© are given 
by the integrals 
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U = 

B = 

h = 
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D x = 
D 3 = 
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dr (g^L^ + ^L? J dr'u d |$ 



'J 2 
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where 



HoH 2 1 — 3 cos 2 
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(10) 

(11) 

(12) 
(13) 
(14) 
(15) 
(16) 

(17) 



is the microscopic dipole-dipolc interaction and = 
$ L( - fi )(r'). Note that the only parameters that appear 



in eqs. ([5]) and ((5]) are A, a and e, which depend on 
both dipolar and s-wave interactions. To compute the 
integrals (fTOj) — (fT6 |) and therefore to determine the param- 
eters (J?])-©, we have used the modes defined as [28j : 
®L(R)(r) = ($ s (r) ± $ as (r))/V2, where $ 8 (r) and $ as (r) 
are the symmetric (ground state) and antisymmetric (first 
excited state) wave functions of the double well potential 
of fig.[T] For the case analyzed here, A = 17.84, a — 0.1307 
and e = 0.1156. 

In the framework of the two-mode model it is possible 
to find an analytical expression for the critical population 
imbalance that gives the onset of self-trapping (see, for 
instance, refs. [TTJ[28]). For 4>{0) = 0, it reads: 



Z P . = 2 



y/1 - aVA" 



A-e 



(18) 



which yields Z c = 0.43 in the present configuration. 

We have obtained the dynamics within the TMM by 
solving cqs. §5§ and ©. The results are given as dotted 
lines in figs. [5] and [31 for the same initial conditions as the 
corresponding TDGP calculation. We see from the figures 
that the TMM is a good qualitative approximation to the 
full dynamics given by the TDGP eq. ([2]): the order of 
magnitude of the frequency and amplitude of the oscil- 
lations, as well as the dynamical regime imposed by the 
initial conditions, arc well predicted. 

The discrepancy between the TDGP dynamics and the 
TMM can be attributed to a combination of different fac- 
tors. Firstly, it has been recently shown that the fact 
that the system does not lie in the deep weak-link limit 
gives rise to a frequency in the two-mode approximation 
larger than the experimental one [3D], which lies closer 
to the TDGP result. Secondly, the self-induced nature of 
the double well means that it depends on time, which is 
not taken into account in a TMM with time-independent 
parameters. Lastly, the SIJJ presented here is clearly 
two-dimensional, whereas the TMM mimics it as one- 
dimensional. Dynamics in other directions different than 
x might affect the behavior of the imbalance and the phase 
difference [291I3T] . This opens the possibility of new phys- 
ical aspects of a SIJJ which are not present in a usual 
J J. 

Extension to other systems. — The whole sys- 
tem is scalable in terms of the dimensionless constants 
A = u! z /u!±, (To = <To/a±, C = AnNa/a±, D = 
Nfj,o/j, 2 m/('l'Kh' 2 a±) and Vo = Vo/{htJ±), which are the co- 
efficients of the different terms of the dimensionless TDGP 
equation. For the case considered here A = 11, cto = 2.08, 
C = 96.77, D = 24.99 and V = 30. With these dimen- 
sionless constants it would be easy to export the same 
physics to another set of parameters experimentally ac- 
cessible in 52 Cr pT2]113|J 32] . or even to condensates of alkali 
gases such as 39 K or 7 Li, where dipolar effects have been 
observed in Refs. [33J and [3J, respectively, making use of 
a soft zero crossing of the scattering length near a Fesh- 
bach resonance. To give an example of an atomic species 
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with a smaller dipolar moment, we consider a condensate REFERENCES 
of 39 K (fj, ~ 1 hb) with 10 6 atoms. The self-induced dou- 
ble well configuration could be obtained in a toroidal trap 
with lj± — 65 x 2ir s _1 , uj z = 715 x 2ir s" 1 and cr = 4/zm, 
which can be achieved experimentally, and with a scat- 
tering length of a = 0.3 ag, which is compatible with 
the 0.06 as resolution reported in rcf. [33] . The period 
of Josephson oscillations would be in this case close to 
100 ms which is well within the lifetime of the conden- 
sate. Moreover, this scenario can possibly be extended to 
other systems, such as exciton-polaritons in semiconduc- 
tors, which can be polarized and where the ring-shaped 
geometry is experimentally easier to implement than an 
external double well [35] . 

It is worth stressing that throughout this letter the two 
junctions have been considered to behave in phase. This 
is the reason why the whole system can be thought of 
as a single SIJJ and analyzed in terms of a two-mode 
model. However, there exists another regime where the 
two junctions behave with opposite phases (out-of-phase 
situation). It corresponds to atoms tunneling through the 
two junctions with opposite velocities. In a more general 
situation the two phase differences would be independent, 
opening thus the possibility of establishing the bosonic 
analog of the superconducting quantum interference de- 
vice (SQUID) pp. This analysis is however beyond the 
scope of the present article. 

Conclusions. — In conclusion, we propose a dipolar 
self-induced bosonic Josephson junction. This junction is 
created by the anisotropic character of the dipolar inter- 
action modulated by a toroidal trap, which gives rise to 
a ring-shaped, double-well effective potential. The time- 
dependent Gross-Pitacvskii equation predicts Josephson 
oscillations as well as a self-trapping regime in this system, 
depending on the initial population imbalance. Moreover, 
the self-induced Josephson junction can be analyzed in a 
two-mode picture, giving qualitative agreement with the 
Gross-Pitaevskii results. The physics of the self-induced 
bosonic Josephson junction is not restricted to chromium 
condensates, but can also be exported to condensates of 
alkali atoms with controllable contact interactions as well 
as to other systems such as exciton-polaritons. We believe 
that it can also be a starting point to study a bosonic 
SQUID interferometer. 
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